Statistics and scaling properties of temperature field in 
symmetrical non-Oberbeck-Boussinesq turbulent convection 

Yuri Burnishev and Victor Steinberg 
Department of Physics of Complex Systems, 
The Weizmann Institute of Science, Rehovot 76100, Israel 
(Dated: September 1, 2011) 

Abstract 

The influence of symmetrical non-Oberbeck-Boussinesq (SNOB) effect on statistical and scaling 
properties of temperature field in turbulent convection is investigated experimentally in SFq in the 
vicinity of its gas-liquid critical point (CP). The main conclusion of the studies is that besides the 
strong Ra and Pr dependence of the rms of temperature fluctuations normalized by the temperature 
difference across the cell, different from the Oberbeck-Boussinesq (OB) case of turbulent convection, 
all rest of statistical and scaling properties of temperature field discussed in details are the same as 
in the OB case. 

PACS numbers: 47.27.-i, 44.25. +f, 47.27.Te, 47.20.Bp 
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I. INTRODUCTION 

Our recent publication {l| introduces a new type of Rayleigh-Benard turbulent convection, 
namely symmetric non-Oberbeck-Boussinesq (SNOB) turbulent convection. The latter is 
characterized by strong temperature and density dependencies of thermodynamic and kinetic 
properties of a supercritical fluid near its gas-liquid critical point (CP) at the average critical 
density of the fluid in the cell. This strong temperature and density dependence results in 
a strong but symmetric height dependence of the main physical properties, which enter 
into the expressions for the control parameters of turbulent convection. So, in spite of 
strong variations of the fluid properties across the cell height, up-down symmetry of the 
temperature drops across the top and bottom halves of the cell and of the top and bottom 
thermal boundary layer widths is preserved. Thus, it was shown [l| that in this case the same 
scaling of Nu = F(Ra, Pr) with the Rayleigh number, Ra = (3gL 3 A/b>K, as in the Oberbeck- 
Boussinesq (OB) turbulent convection is preserved but a much stronger the Prandtl number, 
Pr = v/k, dependence of Nu is found. Here g is the gravity acceleration, (3 is the fluid 
isobaric thermal expansion coefficient, v and k are the fluid kinematic viscosity and thermal 
diffusivity, respectively, Nu = QL/XA describes the heat transport by turbulent convection, 
Q is the heat flux density, L is the the cell height, A is the thermal conductivity of the fluid, 
A = 7], — T t is the temperature difference across the cell, and and T t are the temperatures 
of the bottom and top plates, respectively. In order to single out the influence of the non-OB 
effect on the heat transport, the heat transport experiments we repeated for each Pr with 
an eight-fold larger non-OB effect by decreasing the cell height at the same Ra. As the 
result, the Nu = F(Ra, Pr) scaling was not altered. Therefore, the conclusion has been 
made that the strong SNOB effect by itself is not responsible for the strong Pr dependence 
of the heat transport near CP observed in the experiment but, probably, strongly enhanced 
compressibility, which accompanies an increase of Pr due to approach to CP, could cause 
the observed Pr dependence. 

Next natural question arises whether the SNOB effect results in the same scaling relations 
as in the OB case for statistics of temperature and velocity fields. In our early papers 



3j, we have already studied the statistics, frequency power spectra and scaling of velocity 
and temperature fields. However, this old experiment was less controlled and had lower 
resolution in temperature measurements with a narrower range of Pr and Ra than in the 
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current experiment. The goal of the experiments presented in this paper is to find out 
whether the SNOB effect in turbulent convection modifies the scaling properties of the 
characteristic frequencies, corresponding to large scale circulation, to the Bolgiano length 
and to the dissipation scale, of the temperature power spectra and correlation functions, and 
the scaling properties of the structure functions of the temperature increments compared 
with the OB case. These statistical properties of the temperature field are studied in a wide 
range of Pr and Ra, the same as was reported in our recent paper on the heat transport jl|. 

The outline of the paper is as follows: In Section II the experimental setup is described 
in details. In Section III the local temperature measurements and statistics of temperature 
filed are presented. Frequency spectra of temperature fluctuations and auto- and cross- 
correlation functions of temperature field are discussed in Section IV, and temperature 
structure functions at high values of Ra and Pr are shown in Section V. The results are 
summarized and discussed in section VI. 



II. EXPERIMENTAL SETUP 



As convecting fluid, high purity (99.998%) gaseous sulfurhexafluoride SFq was used in 
the single-phase region in the vicinity of the gas-liquid critical temperature T c and at the 
critical density p c . This fluid was chosen due to its relatively low T c and critical pressure P c , 
and its well-known thermodynamic and kinetic properties, both far away and in the close 
vicinity of the CP. The region of the phase diagram exploited in our experiment is located 
at Pav — Pc and in the region of the reduced mid-plane temperature, (T m — T c )/T c , between 
1.6 x 1CT 4 and 3.4 x 10~ 2 . As we pointed out in Ref. fl, NIST Q gives the following values 
for the critical parameters: T. = 318733 K, P. = 37.5455 bars, * = 743.81 which we 

further use in our analysis. In Ref. we provided detail information about thermodynamic 
and kinetic properties of SFq in the region of our interest and the way of temperature and 
pressure calibration and determination of the critical density p c . 

The experiments were carried out in two cells: one cell of a square cross section a x a = 
70 x 70 mm 2 and height L = 90 mm (aspect ratio T = L/a = 1.286), and one cylindrical 
cell with diameter D = 90 mm and height L = 90 mm (r = L/D = 1). The side walls of 
all cells were made of 4 mm thick plexiglass (X p i = 0.2 W/m K). The cells were enclosed 
into a high-pressure stainless steel vessel. The interior of the cell communicated with the 
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gap between cell and pressure vessel, through a small hole of about 1 mm located in the cell 
mid-plane. For the cylindrical cell, this gap was narrower than 5 mm and tightly filled with 
glass wool (X gw = 0.04 W/m K), tightened by a band. The cell was then inserted into the 
vessel with a sufficient force. In the square cross section cells, the outer space was larger, but 
also tightly packed with glass wool. In both cases the filler was used to reduce the amount 
of gas surrounding the cell, in order to prevent convection outside the cell and to reduce the 
heat transfer due to the outer gas (see Fig. 1). 

In the rectangular cell, a top plate made of h t = 20 mm thick stainless steel 304 
(X ss = 25.8 W/m K) and bottom plate made of = 8 mm thick aluminium 7075 
(Xai = 133 W/m K) were used. In the cylindrical cell, copper plates for both top and 
bottom were used. Heat was supplied at the cell bottom by a metal-film heater of resistance 
100 Q, that covered uniformly the entire active area of the bottom plate. The experiment 
was carried out at a constant heat flux, which was measured during the experiment. 

The high-pressure vessel itself was immersed into a water bath. The top plate was cooled 
down by water circulating in a two loop temperature controlled system. The first water loop 
was refrigerated by a circulator (Lauda Inc.) with a temperature stability of ±10 mK. The 
second water circuit, coupled via a heat exchanger to the first one, cooled down the pressure 
vessel. The water temperature in the second loop was controlled by a heater located just 
before the inlet to the apparatus, computer-controlled by a feedback loop using a thermistor 
placed in the mixing chamber of the water bath. The entire apparatus was covered with 
insulating polyuretane sheets to reduce thermal losses. The temperature stability of the 
bath achieved in a such way was better than 0.3 mK rms. All heat transfer and temperature 
measurements were carried out at fixed volume of the gas at the critical density. To perform 
the local temperature measurements, first the average temperature of the water bath, 
was adjusted above T c . The bottom temperature in the cell, was then increased and 
Tbath decreased in a such way that the temperature in the cell mid-plane T m = (TJ, + T t )/2 
remained constant. This way, the value of Pr at the cell midplane also remained constant 
with an accuracy better than 1% for a given temperature difference A = — T t across the 
cell or for a given Ra. Each time, the value of T t was recalculated from and compared 
with the temperature measured by the thermistor glued on the surface of the top plate. 

Temperature measurements inside the cell were conducted using several thermistors and 
one platinum thermometer of 1 kVt resistance. The latter provided the thermodynamic 
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scale which defines the critical parameters of the fluid, entering in either the parametric 
equation of state \^ or the tabulated data from NIST [5|. Two glass-embedded stable 
thermistors P20BB204M from Thermometries Inc. were used in the water bath and inside 
the bottom plate. Three glass bead thermistors B05KB204M of 150 /im diameter from 
Thermometries Inc. were used to measure temperature inside the cell: one was glued close 
to the center of the top plate and two were mounted on the vertical motor-driven probe at 
the cell center. All thermistors together with the resistance thermometer were calibrated in 
the thermally- regulated water bath (with a temperature stability of ±10 mK) against the 
secondary standard platinum-resistance thermometer of 1 kfl from the Russian Institute 
of Standards. This thermometer was supplied with a calibration accurate to ±10 mK, on 
the ITS90 temperature scale. Long-term stability (during 10 hours) of the thermistors was 
better than 0.1 mK. 

Pressure was measured with a calibrated pressure transducer TJE/727-23 (0-1000 psi) 
from Sensotec Inc. with a full scale output of 4.9975 V, which was read by a 6.5 digits 
multimeter. The pressure transducer was calibrated against the absolute high precision 
pressure gauge Heise Inc, USA (100.00 bar) with accuracy ±25 mbar in the whole range. 
Long-term stability (during 10 hours) of the pressure transducer was better than 0.1 mbar. 

The determination of the critical density of SFq, which is crucial for the experiment near 
CP on the critical isochore, was based on the equation of state J^l and on the NIST data 
We adjusted the gas density far away from CP but in the region still covered by the equation 
of state and the NIST data [4, 5|, by variation of T and P and their precise measurements to 
be equal to p c according to either the equation of state or the NIST data Jj]]. This method 
has an error in p c due to the errors in T and P of the order of 0.4%, while the accuracy in 
determination of T c and P c via the equation of state was completely defined by the errors in 
the absolute values of T and P. 



III. LOCAL TEMPERATURE MEASUREMENTS AND STATISTICS OF TEM- 
PERATURE FIELD IN RAYLEIGH-BENARD TURBULENT CONVECTION 

Local temperature measurements inside the cell were conducted by two glass-bead ther- 
mistors B05KB204M of 150 /im diameter and fast response time of 5 msec (immersed in 
water) from Thermometries Inc., mounted on a stainless steel capillary tube of 0.8 mm outer 
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diameter, separated vertically by a fixed distance of about half of the cell height (about 45 
mm). Each thermistor is mounted on the copper leads of 0.1 mm diameter and 5 mm 
long and suspended on its contacting leads of 0.02 mm diameter and about 5 mm long, 
made of platinum alloy (A = 71.6 W/m K) (Fig. 1). The probe was driven vertically by 
a computer-controlled stepping motor at the cell center in both cylindrical and rectangular 
cells, allowing precise computer control of the thermistor positions. Besides in the rectan- 
gular cell another probe was used at one of the cell corners. Temperature of the thermistors 
on the probe was measured in two ways. One way was the resistance measurements via an 
ac bridge and Lock-in amplifier (PAR model 124) at 83 Hz and A-to-D card. The frequency 
was chosen to be about twice larger than the expected highest frequency in the temperature 
frequency spectra at the highest values of Ra ~ 10 16 , that is about 30-40 Hz. Another way 
was the thermistor resistance measurements via a 6.5 digits multimeter Keithley 2010 with 
a sampling rate of 54 Hz. The comparison showed a good agreement between two data sets 
measured by two techniques. At each value of Ra and Pr the temperature measurements 
were taken at up to 30 locations but the data only at three heights are presented in details: 
10, 22.5 and 45 (center) mm with large statistics in average up to 10 7 data points for each 
set. In the cylindrical cell, the data for 6 different Pr = 11, 25, 74, 130, 180, 244 and up to 5 
values of Ra for each Pr value were taken. So all together 66 data sets for different values of 
Pr, Ra, and 3 locations (10, 22.5 and 45 mm) for each set of Pr, Ra were analyzed. In the 
rectangular cell, the data only for 3 values of Pr = 4.8, 12, 20 and 5 values of Ra for each Pr 
and for the same three locations as above were analyzed. Figure 2 presents segments of tem- 
perature fluctuations time traces ST(t) = T(t) — T for different Ra and Pr at a mid-height 
and at the center of the cylindrical and rectangular cells. Here T is the time average temper- 
ature. It can be seen that positive and negative spikes are distributed approximately equally. 
It is naturally reflected in almost symmetrical probability distribution functions (PDF) of 
temperature fluctuations at the central location in the both geometry cells with exponential 
tails taken at the same values of Ra and Pr as in Fig. 2 (Fig. 3). To compare PDFs at 
different heights at the cell center, we present PDFs of temperature fluctuations at several 
values of Pr and selected values of Ra in the both cylindrical and rectangular cells in Fig. 
4. As one expects, out of the mid-height PDFs show significant skewness and exponential 
tails, as a rule, in the direction of positive temperature fluctuations. The rms temperature 
fluctuations normalized by the temperature difference across the cell, T rms /A, and measured 
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at the cell mid-height as a function of Ra at different values of Pr are shown in Fig. 5. The 
data are fitted by the power law T rms /A = (46 ± ll)Ra a Pr 13 , where a = —0.43 ± 0.03 and 
j3 = 1.24±0.01. The Ra dependence is much stronger than a ~ —0.145 found in the OB case 
6|-|8| while the Pr dependence was never reported. The third M 3 and fourth M4 moments 
of the temperature distributions at the cell mid-height as a function of Ra for different Pr 
are presented in Fig. 6. Significant scatter between zero and unity is observed at M3 and 
between 3 and 4 for M4 without visible dependence on either Ra or Pr that witnesses on 
systematic deviations from the Gaussian and almost symmetrical distribution. Figure 7a,b,c 
shows a dependence of three moments of PDFs of temperature fluctuations on height h at 
Pr = 130 and three values of Ra: 4.6 x 10 13 , 1.9 x 10 14 ,6.1 x 10 14 . The normalized rms 
T rms /A is independent of h at Pr = 130 and three values of Ra as presented in Fig. 7a. 
The same conclusion can be made for all values of Ra and Pr explored in the experiment. 
The skewness (third moment) M 3 and flatness (fourth moment) M 4 show positive deviations 
towards 1 and 4, respectively, at the top and bottom of the cell, as one can see in Fig. 7b, c 
for Pr = 130 and three values of Ra. The same conclusion can be reached for all values of 
Ra and Pr explored in the experiment. 

IV. TEMPERATURE POWER SPECTRA, AUTO- AND CROSS-CORRELATION 
FUNCTIONS OF THE TEMPERATURE FIELD 

Further by using frequency power spectra and auto- and cross-correlation functions of 
temperature fluctuations we are going to study scaling of the frequency of the large scale 
temperature oscillations with Ra and Pr. These frequency peaks in the power spectra 
and the corresponding characteristic oscillation periods in the correlation functions were 

nnn 

reported in many studies on OB turbulent convection |6|, l9Ml 1| . We concentrate first on the 
temperature measurements at the center and the mid-height of the cell in a wide range of Ra 
and Pr. As an example, we present in Fig. 8 the frequency power spectra of temperature 
P(f) at Pr = 126 and 3 values of Ra = 3.1 x 10 13 , 1.1 x 10 14 ,4 x 10 14 and at h = 45 mm in 
the cylindrical cell. A striking feature of all frequency power spectra studied is an emergence 
of either one (in a rectangular cell one observes either one main peak or the main peak and 
a much smaller second harmonic peak) or two sharp peaks at low frequencies. The presence 
of the second harmonic peak at 2/ c , besides the main peak at f c , indicates a nonlinear 
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nature of the oscillations. The second peak can be comparable and sometime even larger (in 
cylindrical cell) than the main one (see Fig. 8). In our early measurements, the main peak 
in the frequency power spectra for temperature fluctuations was barely observed, probably, 
due to lower resolution and smaller statistics, whereas similar peaks in the frequency power 
spectra of velocity fluctuations were detected and analyzed The normalized main peak 
frequencies for both cells at all values of Pr and Ra at h = 45 mm are scaled as f c L 2 /u = 
(0.53 ± o.05) J Ra a43±aoi Pr- a74±0 03 (see Fig. 9) with the coefficient found from the high 
resolution compensating plot, which is presented in the inset in Fig. 9 (we used the value 
v = 5 x 10 -4 cm 2 /s, see Ref. jl|). Good agreement with the result on scaling of the 
frequency of the oscillations found from the peak in the velocity power spectra reported 
in Ref. |3| is found. Similar scaling relations were obtained for the main peak f c L 2 /u in 
the frequency power spectra at the cell center and other two heights below the cell mid- 
plane h = 22.5 and h = 10 mm in the same range of Pr and Ra. As shown in Fig. 10 
a,b, the corresponding scaling relations are f c L 2 /u = (0.3 ± 0.05)Pa a45±a02 Pr _0 ' 74±0 ' 07 and 
f c L 2 /u = (1 ± 0.15) J Ra a41±a02 Pr" a74±a07 at h = 22.5 and h = 10 mm, respectively, with 
close values of exponents for all three heights. Figure 11 shows an independence of f c L 2 /u 
on the height at Pr = 11, Ra = 3 x 10 11 and Pr = 74, Ra = 1 x 10 14 , as examples. 

Similar information can be found from the temperature auto-correlation function (ACF) 
g a (r) = lim-f^oo Jq 8T(t)5T(t + r)dt/A 2 as a function of the delay time r. Here 8T(t) = 
T(t) — T is the temperature fluctuation and T is the average local temperature. Figure 12 
presents several examples of the measured ACF at h = 45 mm and at the following values 
of (Pr,Ra): (25,2.1 x 10 12 ), (130,5.8 x 10 14 ), and (244,1.1 x 10 15 ). First, the measured 
ACF exhibit decaying oscillations. Second, the characteristic oscillation period r c strongly 
depends on Pr and Ra. And third, the characteristic oscillation period r c is equal to 
the inverse frequency f c of the main peak in the corresponding power spectrum, and the 
exponents in the scaling relations of T~ x L 2 jv and f c L 2 /u with Ra and Pr are the same 
inside the error bars. 

Similar features are also observed in the temperature cross-correlation function (CCF) 
g c (r) = limx_>oo / Q T 5Ti(t)8T 2 (t + r)dt/A 2 as a function of r. In the measurements, the two 
thermistors separated by a fixed distance of 45 mm apart are located on the same vertical 
temperature probe, which is driven vertically at the cell center. Figure 13 presents CCF at 
one location h = 22.5 mm of the lower thermistor (and correspondingly of the upper one at 
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h = 67.5 mm ) and at three values of Pr and Ra: (11, 4 x 10 12 ), (74, 3 x 10 13 ), (180, 2 x 10 15 ). 
Similar to the auto-correlation functions, CCFs exhibit oscillations with the frequency r" 1 
equal to the main peak frequency f c in the corresponding power spectrum. The oscillations 
have a decay with the coherence time r^ ec , which is defined from the fit and is about an 
order of magnitude larger than the corresponding oscillation period r c . Scaling relation 
of the normalized oscillation frequency function of Pr and Ra, found from the 

cross-correlation functions at h — 22.5 mm is presented in Fig. 14. The power-law fit gives 
r ( T 1 L 2 jv — (1 ± 0.08)-Ra a41±a02 Pr~ a72±a02 with the exponent values in a good agreement 
with those found for the main frequency peak in the power spectra (see Fig. 10a). We were 
not able to get the scaling relation for the coherence frequency r^L 2 jv due to large scatter 
and insufficient number of the data points. On the other hand, one finds that the ratio 
Tdec/Tc grows from approximately 3.5 at Pr = 11 and Ra — (1 -j-4) x 10 12 till 4.5 at Pr = 74 
and Ra = (3 x 10 13 1 x 10 14 ) and up to 8.8 at Pr = 180 and Ra = (4 x 10 14 4- 2 x 10 15 ). 
It means that the growing number of the observed oscillation periods during the decay 
period indicates the growing coherence of the oscillation of the large scale circulation. As 
shown in Ref. jll|, the large scale circulation occurs before the threshold for the coherent 
oscillations. There is a threshold for the onset of the coherent oscillations, which are resulted 
from coherent action of the rising and falling plumes jll|. Thus, our measurements show 
that the threshold value of Ra for the onset of the coherent oscillations also depends on Pr. 

Another feature of the cross-correlation function is the time delay (shift) T s h of its peak 
value from zero. Examples of the time delays are shown in Fig. 15 a,b for Pr = 244 and 
three values of Ra and for Pr = 130 and another three values of Ra. This time delay 
is probably related to the large scale circulation, and so can supply information about its 
velocity. The Ra and Pr dependence of the normalized time delay is presented in Fig. 16, 
where the corresponding scaling relation is t~^L 2 jv = (12.1 ± 3.g)/£ a - 34 ± - 02 p r -°- 71±0 - 01 
with the coefficient determined from the high resolution compensating plot (see inset in Fig. 
16). Relatively large scatter of the data, compared for example with the data in Figs. 9, 
10, and 14, which is reflected in large error in the constant, could be probably the reason 
for the smaller value of the exponent in Ra dependence 0.34, compared with that obtained 
for the normalized oscillation frequency t c L 2 /u of the cross-correlation functions 0.41 (see 
Fig. 14 and the corresponding scaling relation above). As shown in Ref. jll|, the large 
scale circulation velocity and the oscillation frequency are related to each other by simple 
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relation, which in frequency domain gives the ratio 4. Due to different scaling exponents for 
Ra dependence mentioned above, the ratio between and r c _1 varies between 1.3 and 2.2, 
which difference from 4 can also be caused by the large scatter and insufficient statistics. 

To further examine the scaling properties of the frequency power spectra of temperature 
fluctuations we use the method introduced for hydrodynamic turbulence and suggested for 
convective turbulence in Ref. 12|. The method utilizes the peak frequency f p corresponding 
to the maximum value of the temperature dissipation spectra f 2 P(f) as the characteristic 
frequency to collapse the power spectra for different Ra and Pr. Figure 17 shows scaled 
frequency power spectra of temperature fluctuations P(f)/P(f p ) versus f / f p for all values of 
Pr and Ra explored in the experiment in the cylindrical cell. An example of the dissipation 
spectrum f 2 P(f) at Pr = 244 and Ra = 1.6 x 10 13 at the height h = 45 mm and the cell 
center is presented in Fig. 18. The solid line in Fig. 18 is the sixth-order polynomial fit 
to the data, which allows an accurate determination of the frequency f p corresponding to 
the data maximum. It is clear from Fig. 17, that for each Pr and Ra there is well-defined 
f p , with respect to which the frequency power spectrum is a universal function. Thus, we 



checked that the suggested method 



12j works well in a wide range of Pr and Ra up to 



their very high values. To compare the scaling results for f c and f p , we present in Fig. 19 
the normalized frequency f p L 2 /u as a function of Pr and Ra. The corresponding scaling 
relation is f p L 2 /is = (162 ± 21) J Ra a27±0 01 Pr- a43±0 02 with the coefficient obtained from 
the compensating plot shown in the inset in Fig. 19. Thus we found experimentally two 
characteristic frequencies, f c and f p in statistics and dynamics of the temperature field and 
studied their scaling relations in regards to the two control parameters of the problem Pr 
and Ra. 

Finally, the scaled frequency power spectra of the temperature fluctuations exhibit rather 
clear power-law region in the low frequency end with the exponent 1.4 equal to the Bolgiano- 
Obukhov (BO) scaling though in the wave number domain that is more convincingly demon- 
strated by the high resolution compensating plot of P(f)/ P(f p )(f / f P ) 1A versus f/f p in the 
inset in Fig. 17. This is the only clearly identified scaling region found in the scaled power 



spectra. On the other hand, it is rather well established |12l . |13j, that in developed tur- 
bulent convection at large values of Ra and sufficiently small Pr, two scaling regions can 
be observed simultaneously. One of them is the BO region, where buoyancy controls the 
dynamics, and another is the inertial, Kolmogorov range, in which the effects of buoyancy 
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are irrelevant [l^. These regions are separated in scale by the Bolgiano length scale l B , and 
at scales I > Ib the BO scaling is observed, while at I < Ib the Kolmogorov scaling is found 
1^ [ . According to 
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15j, Ib can be estimated by 



l B = LNu^iRaPr)- 1 / 4 , (1) 

and can be considered as a characteristic value of Ib- Using the expression of Eq. (1), we 
can estimate the characteristic Bolgiano frequency based on large scale circulation velocity 
defined via f c . Then the Bolgiano frequency can be written in the normalized form as 



f B {L 2 /v) = (L 2 /u)f c L/l B = 0.64Ra a53 Pr-°- 39 , (2) 

where Nu = 0.18Ra°' 3 Pr~ ' 2 , taken from Ref. is substituted into the expression of 
Eq. (1), and the expression for f c is used from the fit in Fig. 9. We point out that in the 
case considered the scaling of fs in both Pr and Ra is very different from the above found 
scaling for f p in contrast to Ref. Jl2j. The ratio Jb/ fp = 3.9 x 10~ 3 i?a ' 24 change rather 
strongly with Ra but practically independent of Pr in the range of Pr and Ra considered 
in the experiment. Finally, the largest characteristic frequency in the temperature power 
spectra is the dissipation one defined from the dissipation (or Kolmogorov) length 14 



l d = LiPrRaNu)- 1 ^ = 1.54Ra- a325 Pr-°- 2 , (3) 

with Nu substituted from Ref. Jl|. The corresponding dissipation frequency is f^is/L 2 ) = 
0A2Ra°' 75 Pr~ 5A . In Fig. 17 we show by arrows for comparison the minimal and maximal 
values of f c j f p and fs / f P for some specified values of Pr and Ra in the whole range of these 
parameters studied, while fd/f p remains outside of frequency range presented in the plot at 
all values of Pr and Ra. From our estimates one finds that /# is located rather close to 
the intersect of P(f) with the noise floor and therefore rather close to the cut-off frequency 
of the frequency power spectra for the most of the values of Pr and Ra considered in the 
experiment. Only for the lowest values of Pr and Ra one can expect to find f B close to 
f p (see Fig. 17). To conclude, the experimental results for f c and f p and the theoretical 
estimates for Jb and fa in the range of Pr and Ra considered in the experiment give the 
following range of changes for these characteristic frequencies: f c — [1.1 x 10~ 2 -j- 0.13] Hz, 
f p = [0.3 -=- 1.5] Hz, f B = [1 -=- 75] Hz, and f d = [45 3 x 10 4 ] Hz. So, f d is out of scale of 



11 



the power spectra obtained in the entire range of Pr and Ra, where the cut-off frequency is 
exclusively determined by the noise floor, and fs for low Pr and low Ra probably could be 
found at frequencies significantly above the intersection with the noise floor. 

To proceed further with the analysis of an experimental determination of and probable 
existence of the mentioned above two scaling regions we use a different approach based on 
temperature structure functions. 



V. TEMPERATURE STRUCTURE FUNCTIONS IN THE BOLGIANO- 
OBUKHOV REGIME OF THERMAL CONVECTION AT HIGH Pr AND Ra. 

According to the values of f c and f p obtained experimentally and the values of fs ob- 
tained from estimates, one can expect at some values of Pr and Ra from the range of the 
values explored in the experiment to find /# rather close to f p , between f c and the cut-off 
frequency, and therefore to observe the second scaling region at f> fn corresponding to 



the Kolmogorov scaling. As was shown in the recent publications |16|. |17|. a generalized ex- 
tended self-similarity (ESS) method based on structure functions of temperature increments 
can be used to analyze the temperature data. Let us first to examine the behavior of the 
structure functions S p (t) = (\T(t + r) — T(t)\ p ) up to order p = 8. Similar to the correlation 
functions discussed in the previous Section, the structure functions at least up to p = 4 
exhibit oscillations (Fig. 20), which minima are associated with one, two, etc. coherent 
oscillation periods of the large scale circulation. Naturally, the oscillation period and the 
scaling with Pr and Ra are the same as for the auto- and cross-correlation functions and 
the inverse frequency f c found from the power spectra (see Fig. 21). 

To proceed further with the ESS method, let us use the approach, suggested in Ref. 1X8 1 
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and employed for the temperature data analysis in turbulent convection in Refs. 
and plot the ratio Si(t)/[S2(t)] 1 ^ 2 versus time increments r, as it is presented in Figs. 22 
and 23 for several values of Pr and Ra. On the plots in Fig. 22 a,b the data at (a) Pr = 11 
and i?a = 4 x 10 11 and (b) Pr = 74 and Ra = 1 x 10 14 are shown with scaling in the entire 
BO region (indicated by the solid line) between r c and tb (found from the estimates) and 
with oscillations on scales r > r c . The corresponding scaling exponent is close to 0.4 on the 
plot S*2(r) (not shown) that corresponds to the exponent 1.4 found from the temperature 
power spectra (see Fig. 17). On the plots in Fig. 23 a,b for the data at (a) Pr = 180 and 
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Ra = 2 x 10 15 and (b) Pr = 244 and Ra = 3 x 10 15 a change in the slope occurs at r B = 0.59 
s and tb = 0.93 s, respectively, the time scales close to these found from the estimates. Two 
scaling ranges are identified above and below the breaking point, which we associate with 
Tb, corresponding to the Bolgiano length lg. Thus, at r > tb the BO scaling regime is 
identified, while at r < Tb the Kolmogorov scaling regime is found. 

Figure 24 shows the structure functions S p up to order p = 8 in the ESS presentation S p 
versus S2 for Pr = 130 and Ra = 2 x 10 14 in the BO range of scales. Reasonable scaling 
is observed though it deteriorates slightly for the higher order structure functions. The 
normalized scaling exponents Cpl C2 obtained from the plots in Fig. 24 are shown in Fig. 
25 in comparison with the results of averaging over all the data for all Pr and Ra studied 



3 



taken on turbulent convection in 



20] and in 



in the experiment. We also show the data from Ref. 

helium, the data on passive temperature fluctuations in turbulent flows in air 
helium 2l| and theoretical suggestion 22[. All results exhibit strong intermittent behavior 
of the temperature field with our data found surprisingly close to the theoretical curve 
corresponding to a passive scalar analysis. 

VI. DISCUSSION AND CONCLUSIONS 

The detail studies of the statistical and scaling properties of the temperature fluctuation 
field of SNOB turbulent convection in a wide range of Pr and Ra show that they are rather 
close to those investigated in the OB case. The only difference we found is in much stronger 
Ra dependence of the normalized rms temperature fluctuations T rms /A than in the OB 
case and unusually strong Pr dependence, which cannot be compared with the OB case 
due to lack of the data. On the other hand, scaling of the peak frequency of the coherent 
oscillations in the temperature power spectra as well as the corresponding periods of auto- 
and cross-correlation functions and structure functions agree well with the those in the OB 
case. Another visible difference found is the emergence of strong second harmonic peak 
sometimes even higher than the main one that points out on the strong nonlinearity of the 
coherent oscillations. We also found that the degree of coherence of the oscillations depends 
on Pr besides the Ra dependence studied in the OB case in Ref. 



the method suggested in Ref. 



111 ]. We verified that 



12| to collapse the temperature power spectra using the the 



peak frequency /„ of the corresponding dissipation spectra works well in a wide range of 



13 



Pr and Ra in SNOB turbulent convection though the corresponding scaling of f„ differs 
significantly from the estimated scaling of fs in contrast to the results of Ref. (l2|. From 
the structure function analysis we establish that for the most of the range of Pr and Ra 
studied in the experiment only the BO scaling region is observed though for the smallest 
small values of Pr and Ra both scaling regions were identified. And finally, the dependence 
of the normalized scaling exponents of the structure functions ( p / (2 on p agrees rather well 
with those for the OB case and is found surprisingly close to the theoretical predictions for 
a passive scalar behavior 22 1. 
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Figure 1: Schematic drawings of a convective cell and of a thermistor probe for local temperature 
measurements. 
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Figure 2: Time series of temperature fluctuations ST(t) = T(t) — T at the cell center and at mid- 
height h = 45 mm for the following values of Pr and Ra: (a) (11, 4 x 10 12 ); (b) (74, 1 x 10 14 ); (c) 
(180, 2 x 10 15 ); (d) (244, 3 x 10 15 ). 
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Figure 3: PDFs of the normalized by the corresponding rms of temperature fluctuations 5T/T 7 
at the cell center and at mid-height h = 45 mm for the same values of Pr and Ra as in Fig. 2. 
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Figure 4: PDFs of the normalized by the corresponding rms of temperature fluctuations 5T/T rms at 
the cell center and at following values of heights h, Pr and Ra: (a) (10 and 22.5 mm, 11, 4.7 x 10 12 ); 
(b) (10 and 67.5 mm, 74, 1.1 x 10 14 ); (c) (67.5 and 80 mm, 180, 2 x 10 15 ); (d) (5 and 22.5 mm, 
244, 3.7 x 10 15 ). 
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Figure 5: Normalized rms of temperature fluctuations function of Ra for different Pr. 

The solid line is the fit to the data T rms /A = (46 ± ii)jfa-o.43±o.03p r i.24±o.oi_ xhe inset shows 
the compensated plot of (T rms / A) Ra 3 p r ~ 124 versus Ra, from which the constant is found. 
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Figure 6: Third (skewness) M3 and fourth (flatness) M4 moments of PDFs of temperature fluctu- 
ations in a whole range of Pr and Ra studied. 
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Figure 7: Dependence of (a) the normalized dispersion T rms /A, (b) skewness M3, and (c) flatness 
M4 of PDFs of temperature fluctuations on h at Pr = 130 and three values of Ra: 4.6 x 10 13 , 1.9 x 
10 14 ,6.1 x 10 14 . 
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Figure 8: Frequency power spectra of temperature fluctuations P(f) at the cylindrical cell center 
and mid-height h = 45 mm for Pr = 130 and three values of Ra: 3.1 x 10 13 , 1.1 x 10 14 ,4 x 10 14 , 
where the second peak 2/ in the power spectra sometime larger than the main one. 
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Figure 9: Normalized frequency of the main peak in the power spectra at all values of Pr and 
Ra under studies at h = 45 mm in both types of cells. The solid line is the fit for all the data 
f c L 2 /u = (0.53 ± 0.05)«a - 43±0 - 01 Pr-°- 74±0 - 03 . 




Figure 10: Normalized frequency of the main peak in the power spectra f c L 2 jv at all values of 
Pr and Ra under studies at (a) h = 22.5 mm and (b) h = 10 mm in both types of cells. The 
solid line are the fits for all the data (a) f c L 2 /u = (0.3 ± 0.06) J Ra a45±a02 Pr-°' 74±0 ' 07 and (b) 
f c L 2 /u = (1 ± 0.15) J Ra°- 41±0 - 02 Pr-°- 74±0 - 07 . 
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Figure 11: Normalized frequency of the main peak in the power spectra f c L 2 fv as a function of 
height h at Pr = 11, Ra = 3 x 10 11 (lower data), Pr = 74, Ra = 1 x 10 14 (upper data). 




10 20 30 40 




10 20 30 40 50 



t, sec 



Figure 12: Several examples of auto-correlation functions of temperature fluctuations measured at 
three locations h = 45, 22.5, 10 mm and the following values of Pr and Ra: (a) 25, 2.1 x 10 12 ; (b) 
130, 5.8 x 10 14 ; (c) 244, 1.1 x 10 15 . 
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Figure 13: Several examples of cross-correlation functions of temperature fluctuations measured 
between two locations at h = 22.5, h = 67.5 mm and the following values of Pr and Ra: (a) 11, 
4 x 10 12 ; (b) 74, 3 x 10 13 ; (c) 180, 2 x 10 15 . The dotted curve is the full fit and the dashed curve 
is the exponential fit to the decay. 




Figure 14: Normalized frequency t~ 1 L 2 /v obtained from the cross-correlation functions as a func- 
tion of Pr and Ra at h = 22.5 mm. The data based on the cross-correlation functions are compared 
with those obtained from the frequency power spectra at h = 22.5 mm and good agreement is found. 
The solid line is the fit t~ x L 2 jv = (1 ± 0.08) J Ra a41±a02 Pr-°' 72±a02 . 
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Figure 15: Examples of the time delays of the peak of the cross-correlation functions at Pr = 244 
and three values of Ra (upper plot) and at Pr = 130 and another three values of Ra (lower plot). 



o 
o 

+ 
o 



10° 



Q. 

> 10 5 



p 10 





©p-lff 










: 5 „ 10 

Q. 


□ □ □ 


.> 


Ra 


^10" 10 13 10 15 






10 12 „ 10 14 

Ra 


10 16 



Figure 16: Normalized frequency T~^l?jv obtained from the cross-correlation functions as a 
function of Pr and Ra at h = 22.5 mm. The solid line is the fit t~^I? jv = (12.1 ± 
3.6)i?a ' 34±0 ' 02 Pr _0 ' 71±0 01 with coefficient taken from the compensating plot in the inset. 
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Figure 17: Scaled temperature power spectra taken in the cylindrical cell at h = 45 mm for all 
values of Pr and Ra explored in the experiment. Arrows indicate minimum and maximum values 
of fc/fp and fs/fp- minimum value of f c /f p ~ 0.07 at Pr = 74 and Ra = 7 x 10 11 ; maximum 
value of fc/fp ~ 0.18 at Pr = 180 and Ra = 2 x 10 15 ; minimum value of fs/ fp ~ 1-66 at Pr = 11 
and Ra = 9 x 10 10 ; maximum value of f B /f p « 20.2 at Pr = 244 and Ra = 3 x 10 15 . The dash 
line has a slope of 1.4. The inset: the compensating plot P(f)/P{fp){f/fp) 1A versus / / f p . 




Figure 18: An example of temperature dissipation spectrum P(f)f 2 versus / taken in the cylindrical 
cell at h = 45 mm at Pr = 244 and Ra = 1.6 x 10 13 is shown. The solid line is the sixth-order 
polynomial fit to determine the peak frequency f p . 
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Figure 19: Normalized frequency of the peak of the dissipation spectra as a function of Ra at 
various Pr. The solid line is the fit f p L 2 /u = (162 ± 21)i2a a27±a01 Pr- a43±a02 . 
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Figure 20: Temperature structure functions S p up to p = 4 for (a) Pr = 74 and Ra = 1 x 10 14 ; (b) 
Pr = 244 and Ra = 3 x 10 15 . 
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Figure 21: Normalized period of oscillations of the structure functions as a function of Pr and Ra 
in both cells. The solid line is the fit t~ 1 L 2 /u = (0.89 ± 0.16)i?a°- 42±0 - 01 Pr- a73±a02 , where the 
coefficient is defined from the compensated plot presented in the inset. 
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Figure 22: The ratio S*i(t) /[^(t)] 1 / 2 versus time increments r taken in the cylindrical cell at 
h = 45 mm for (a) Pr = 11 and i?a = 4 x 10 11 ; (b) Pr = 74 and i?a = 1 X 10 14 . Solid lines show 
the BO scaling region. 
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Figure 23: The ratio Si(t)/[S2(t)] 1 ^ 2 versus time increments r taken in the cylindrical cell at 
h = 45 mm for (a) Pr = 180 and Ra = 2 x 10 15 ; (b) Pr = 244 and Ra = 3 x 10 15 are presented. 
Solid lines are the fits. Arrows in the insets show in (a) tb = 0.59 s and (b) tb = 0.93 s. Solid lines 
show two different scaling regions: the BO region at r > tb and the Kolmogorov region at r < tb- 




Figure 24: Structure functions S p up to order p = 8 versus S2 in the ESS presentation are plotted 
for the data taken in the cylindrical cell at h = 45 mm at Pr = 130 and Ra = 2 x 10 . 
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Figure 25: The normalized scaling exponents Cp/ C2 versus p for the data taken in the cylindrical cell 
at h = 45 mm at Pr = 130 and Ra = 2 x 10 14 (solid triangles) and for the averaged over all Pr and 
Ra (open down-triangles) in comparison with the data on temperature fluctuations obtained in the 
experiments on turbulent convection in helium (open up-triangles) and taken from Ref. 1 1 911 . with 
the data on passive temperature taken in turbulent flows in air (open squares) from Ref. [20j and 



in helium (solid circles) from Ref. 



2l| . and theoretical result for passive temperature fluctuations 



(solid line) from Ref. |22l | are shown. 
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